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We discuss the temperature correction to the Casimir force between nonideal metallic bodies 
which caused disagreement in the literature. A general method to find the troubling term is proposed 
that does not require a direct reference to the Lifshitz formula. The linear in temperature correction 
is shown to survive for nonideal metals. It is important for small separations between bodies tested 
in the recent experiments. 
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The Casimir force [Q has been measured with high precision in recent experiments |2|-|5j between metallized sphere 
and plate. Also there are plans (If} to look for very weak hypothetical forces where the Casimir force is the main 
background. All this makes the precise evaluation of the Casimir force an important problem. Here we will discuss 
a particular problem concerning the temperature dependence of the force between macroscopic bodies made of real 
metals. Discussion of this problem in the literature revealed disagreement between different authors and at the 
moment there are three different results. 

For perfect conductors the temperature correction has been found many years ago |^-|ic|l and it is small for small 
separations between bodies a -C ch/kT or equivalently for low temperature. For a sphere above a plate the leading 
term behaves as ( T/T e f f) 3 , where kT e ff = hc/2a. This result follows from a general expression for the Casimir force 
given by Lifshitz Jll|,|l2| modified for the case of sphere-plate geometry with the proximity force theorem (PFT) jl3) : 
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e {iQn) is the dielectric function of the used material at imaginary frequencies. The prime over the sum sign indicates 
that the first term n = has to be taken with the coefficient 1/2. Throughout this paper we will work with the force 
F(a) between sphere and plate. The force for the plate-plate configuration can be restored as F' (a) /2ttR. 

For small temperature the sum in (Q) can be replaced by the integral and the resulting force does not depend on 
the temperature at all. This limit has been considered (IJJl^l to calculate the force in the conditions of experiments. 
In general, the replacement is true with the precision ~ T/T e ff. For the atomic force microscope (AFM) experiments 
HQ the smallest separation was 0.1 \im and the replacement error can be as large as 3% at T = 300° K. It exceeds 
the experimental errors ~1% and, therefore, the finite temperature effect has to be taken into account. We define the 
temperature correction AtF as the difference between forces written as the sum over n and as the integral instead of 
this sum. 

The physical discussion of the first term in the sum (|l|) raised heated debate in the literature. It has a long history. 
For ideal metals the temperature correction found from the Lifshitz formula fl6f did not agree with that found with 
different methods Schwinger et al. |To| clearly demonstrated that when the static limit was correctly taken in 

the Lifshitz theory the discrepancy disappeared. The formulated condition is known as the Schwinger prescription: e 
should approach infinity before the frequency is allowed to go to zero. 

Different results have been found for the n = term in the recent works where the Casimir force was calculated 
between nonideal metals. All the results one can parametrize JfTJ with the following relation 
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F n=0 (a) = a—((3), (3) 

where £ (m) is the zeta-function and a is the parameter. In our paper fL8| a = 1 was found using the Drude model 
for e (i and the Schwinger prescription. Doubt in applicability of the prescription for real metals was expressed 
in p9[| where straightforward calculation gave a = 1/2. Using the plasma model for the metal dielectric function, 
a = 1 — 2c/aoj p was found [po|] (see also pi[| ), where uj p is the plasma frequency. Physical arguments based on the 
finite size of plates in real experiments have led Lamoreaux to the result a = 1 |l7| . There is no disagreement for the 
other terms in the sum (]l|), so different a are responsible for different temperature corrections. In this letter we will 
try to clarify the problem. 

The n = term (|J) gives the classical contribution to the force ( ^-independent). At large separations between 
bodies (T e ff <C T) this is the only term in the sum (g) which survives. All the other terms are exponentially 
suppressed. In this sense one can say that (||) describes the contribution of the long wavelength or low frequencies 
fluctuations (the static limit). Below we demonstrate that the value of a can be found for real metals from very general 
analysis without direct reference to the Lifshitz formula (|l|) if one knows its value for ideal metals. It is important 
because the static limit calculated from (0) is controversial: it gives different values for a if the metal is described by 
the plasma or the Drude models and there is no way to reconcile the values in the limit of zero relaxation frequency. 

The force in the classical limit @) with an arbitrary parameter a can be written using only the dimensional analysis. 
More exactly one can reproduce the force between two plates and then apply PFT. The value of a can depend on the 
material parameters which define the metal dielectric function e (i£) at low frequencies. In this limit e (iC/) depends on 
the only one parameter which is the metal conductivity a or equivalently the resistivity p = l/a. The dimensionless 
parameter a can be a function of the only one dimensionless variable 
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where eo is the free space permittivity. For ideal metals p — > and a (0) = 1 as was established many years ago || [To| . 
Let us estimate the length scale L = cpeo for real good metals. Using the typical resistivity p < 10 ptt ■ cm one finds 
L < 3 A. Of course, one can expand the function a (L/a) in the powers of L/a but there is no too much sense in this. 
The reason is that L is a microscopic scale. Really, it is closely connected with the screening length L s ~ vp/to v of 
degenerated electron gas, where vp is the Fermi velocity of electrons. If 5 = c/uj p is the field penetration depth and I 
is the mean free path for electrons, then L ~ L s (5/1). For good metals S ~ I and L is the same order as L s . Taking 
into account such a short scale will be excessive for the macroscopic Casimir force. Therefore, we can conclude that 
the force in the static limit is the same as for ideal metals, i.e. a = 1. Intuitively it seems quite obvious. In the static 
limit there is no difference between metals because the free charges manage to follow the slowly varying field. That 
is why the boundary conditions in this case do not include any characteristic of a particular metal. This conclusion 
means that if we try to calculate the force using the Lifshitz formula ([!]), the Schwinger prescription has to be applied 
to the n — term independently if it is a perfect or real metal. Lamoreaux |l7j arrived at the same conclusion from 
quite different physical consideration. 

In rcf. |l9| fair calculation of the n = term in ([[]) gave a = 1/2. This value does not depend on any material 
parameter and, therefore, should be the same as for the ideal metal. However, we know that for a perfect conductor 
a = 1 due to the Schwinger prescription. So the disagreement has the same nature as in the old story described 
above. 

In ref. [^0| the plasma model was used for the metal dielectric function 

e(0 = l + ^f. (5) 

This model gives good description of a metal in the infrared range but cannot be used at low frequencies where 
e = 1 + (eopC) -1 - Therefore, it will be good everywhere except of the n — term for which one has to use the 
Schwinger prescription. It is no wonder then that in the static limit a significant correction to a — 1 has been found 
p0| : a = 1 — 2c/auj p . That is because in the plasma model the length scale L — c/uj p (penetration depth) becomes 
macroscopic L > 150 A and expansion in L/a has sense. 

The plasma model cannot be used at low frequencies even in principle. The only way to get s (i£) is the dispersion 
relation 



oo 

£ «) = 1 + ^ J duj 



me" H 

LO 2 + C 2 ' 



(6) 



2 



where e" (ui) is the imaginary part of the dielectric function at real frequencies. In the plasma model there is no 
dissipation and e" (u) — 0. Eq. (^J) one can derive from (||) only calculating the integral with some small but finite 
relaxation frequency uj t . For finite dissipation always there is such a frequency for which e (iC) = 1 + (sopC)^ 1 with 
an arbitrary small p. 

The authors pQ ] claim that some identical transformation of the Lifshitz formula removes the controversy in the 
force for the plasma and Drude models without the use of Schwinger prescription. This statement is definitely wrong 
because any identical transformation cannot change the result. Below it will be demonstrated explicitly. 

Let us discuss now the temperature correction to the Casimir force. The sum for n > in (|lj) as a function of 
temperature contains a piece linear in T which exactly cancels for ideal metals the n = term giving in the small 
temperature limit T/T e fj < 1 the well known result || 
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where -Fb( a ) = T^TicRj (360a 3 ) is the bare Casimir force between sphere and plate. 

If we are using the dielectric function of a real metal, the cancellation of the first term in (llj) is incomplete and 



with a = 1 was used for numerical 
the temperature correction at the 



the linear in T contribution survives. That was noted first in [ JL8J , where eq. 
evaluation of the Casimir force. It was found that for the AFM experiments 
smallest separation is 4 pN against the experimental errors 1 pN. We can give now an explicit expression for the 
correction. 

In ref. [po[ the authors did not use the Schwinger prescription in the n = term and took the plasma model for 
e(i£). For this case they convincingly demonstrated that there was no the linear in T correction and the leading 
correction is only (T/T e ff) . As was shown above the prescription must be applied to the n = term so the only 
difference in evaluation of the sum comes from this term. This difference is just the correction we are looking for 
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The integral here can be interpreted as the linear in T term contained in the sum ([!]) for n > and, of course, it 
can depend on the material parameters because the summation is going over nonzero frequencies £ n . On the other 
hand, since this integral appeared as the n = term in (Q), we should take the functions G1.2 at x n = 0. In this 
limit G2 = 1 but G\ ^ 1. Using then the relation J Q drain (I — e~ x ) = —( (3) one finds the final expression for the 
correction linear in T: 
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Let us stress that (^) is true only for the plasma model. When lu p — > 00 the correction disappears as it should be. 
Expansion in powers of (3 gives 

AtF = (3) ' 8/3 (1 " 3/3 + ° ^ + ° (( T / Te ^ 3 ) ■ ( 10 ) 

For uj p — 2 • I0 16 s _1 , a = 0.1 /im and R = 100 fim as in the AFM experiments ||f|] one gets AtF « 2.5 pN using 

tor calculating directly with the help of (|f]) and 2.9 pN using (|f(]). In the condition of the Lamoreaux experiment 
(uj p — 1.4 ■ 10 16 s -1 , a = 0.6 /ira and R = 12.5 cm) the linear in T part of the correction AtF w 29 pN is within 
the experimental errors. This conclusion is in contrast with that of Bostrom and Sernelius p9| due to obvious reason. 
If a = 1/2 as in ref. |li, the n = term is not canceled by the sum of those at n > and there is no an additional 
suppression ~ (3 as in (|10|). 

The correction increases further if we will use the Drude dielectric function 
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In this case it has to be evaluated numerically using ([!]) and subtracting the same formula with the integral instead 
of the sum. The relaxation frequency uj t influences mostly on the integral since it changes low frequency behavior of 
the integrand. For typical value uj t = 5 • 10 13 s _1 we found AtF « 4.0 pN in condition of the AFM experiments 
HQ . This value coincide with that reported before |Q . In the experiment [Q the smallest separation was a = 63 nm 
and the absolute value of the correction is larger AtF w 10 pN [^2| though its relative value decreases. It has to 
be compared with the experimental error 3.5 pN. Note that this correction make the agreement between theory and 
experiment better p2|. 



In ref. 1 20 1 the authors claim that the Drude model has a principal drawback which is discontinuity of G\ at x = 
for the n = term. Namely, G\ = 1 at x = but G\ — at x ^ in the limit xq — > 0. In this connection we would 
like to stress that the finite discontinuity in one point has no effect on the integral in ([!]). As was explained above 
the problem with the n = term has the physical nature. One has to accurately takes into account the static limit 
for which e — > oo must be taken before allowing £ — > for both real and ideal metals. 

The problem with the n = term was declared [0 to be connected with the discontinuity of G\. This has led the 
authors p(J to the wrong conclusion that the problem can be solved by the mathematical methods. They transformed 
eq. (Jl|) with the Poisson formula and represented it as the Fourier transformation 
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It was stated that in this form the Lifshitz formula had no problem for the Drude model. However, the Poisson 
formula is just an identical transformation for the series and it cannot give any new result. As we know without 
the Schwinger prescription it follows from (|l|) that a = 1/2 jljj]. The same value one can get from (p"2|). It was 
demonstrated numerically ]23| and easy to do analytically. 

To see it explicitly, let us consider the static limit t>1. Calculating the sum in @ for finite but large M » t 
one finds 
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The integral is saturated at ttz/t ~ l/(2Af + 1) <C 1 and we can replace sin7rz/r rs ttz/t. Then can be 
represented as 
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We get the expected result from (Oh taking the limit M — ► oo: 
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It coincide with the n = term in (|l|) as it should be. Therefore, using the Drude model we will get a = 1/2 from 
( |l2| ) so as from the original Lifshitz formula ([!]). 

In conclusion, we have considered the linear in temperature correction to the Casimir force at low temperatures or 
equivalently at small separations. Special care has to be taken to get the contribution of the fluctuations in the static 
limit (n — term). This contribution is canceled for ideal mirrors but cancellation is incomplete for real metals. A 
general method based on the dimensional analysis has been proposed to calculate the troubling n — term without 
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direct reference to the Lifshithz formula. On this bases it was shown that the Schwinger prescription had to be applied 
in the static limit independently if it is real or ideal metal. The static contribution is the same for the plasma and 
Drude models. The temperature correction is shown to be important for the AFM experiments |3p|| at the smallest 
separations but is less than the experimental errors for the torsion pendulum experiment 0. 
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